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There has been an ongoing debate about the role of linear instability waves in the prediction of jet noise. 

Parallel mean flow models, such as the one proposed by Lilley, usually neglect these waves because they cause the 
solution to become infinite. The resulting solution is then non-causal and can, therefore, be quite different from the 
true causal solution for the chaotic flows being considered here. The present paper solves the relevant acoustic 
equations for a non-parallel mean flow by using a vector Green’s function approach and assuming the mean flow to 
be weakly non-parallel, i.e., assuming the spread rate to be small. It demonstrates that linear instability waves must 
be accounted for in order to construct a proper causal solution to the jet noise problem. . Recent experimental results 
(e.g., see Tam, Golebiowski, and Seiner, 1996) show that the small angle spectra radiated by supersonic jets are quite 
different from those radiated at larger angles (say, at 90°) and even exhibit dissimilar frequency scalings (i.e., they 
scale with Helmholtz number as opposed to Strouhal number). The present solution is (among other things )able to 
explain this rather puzzling experimental result. 


1. Introduction 

Lighthill (1952, 1954) provided a systematic basis for predicting jet noise when he rearranged the Navier- 
Stokes equations into the form of a linear wave equation for a medium at rest with a quadrupole-type source term 
(which includes a pressure/density contribution that Lilley (1974) showed to be more appropriately described by a 
dipole-type source.) The crucial step in this so-called acoustic analogy approach amounts to assuming that the 
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source term is in some sense known or that it can at least be modeled in an approximate fashion. Early efforts to 
improve the Lighthill approach focused on accounting for mean flow interaction effects. Phillips (1960), Lilley 
(1974), and many others, sought to accomplish this by rearranging the Navier-Stokes equations into the form of an 
inhomogeneous convective, or moving medium, wave equation rather than the inhomogeneous stationary medium 
wave equation originally proposed by Lighthill. Current industrial noise prediction methods, such as GE’s MGB 
approach (Balsa et al. 1978), are based on a form of the convective wave equation proposed by Lilley (1974) in 
which the wave operator is appropriate to sound propagation on a parallel mean flow and therefore possesses 
homogeneous solutions corresponding to spatially growing instability waves on that flow (Betchov and Criminale, 
1967). 

The complete solution to this equation consists of a particular solution plus these homogeneous contributions, 
but the result cannot be used to calculate the far-field noise because the instability waves become unbounded 
(infinite) far downstream in the flow . The usual resolution to this dilemma is to completely neglect the contribution 
of the instability waves. Unfortunately, the resulting solution turns out to be non-causal, which may be particularly 
serious in flows that support instabilities (i.e., chaotic flows) because small changes in initial (and/or boundary) 
conditions can produce large (i.e., 0(1)) changes in the steady state solution. Arguments against imposing a 
causality requirement (Mani, 1976; Dowling, et. al. 1978 ) usually amount to asserting that it is unnecessary 
because it is not possible to identify an initial time before which the fluctuations (about the base flow) have been 
“switched on’ in an acoustic analogy approach and that a boundedness requirement can, therefore, be imposed on 
the solution. 

A better approach might be to begin with an equation appropriate to sound propagation on a non-parallel flow, 
say the actual mean flow in the jet. The most important difference between this approach and Lilley’s parallel flow 
result is that the homogeneous solutions to the acoustic equations correspond to instability waves that grow and then 
decay on the diverging, non-parallel base flow and therefore always remain bounded — which eliminates the 
dilemma alluded to above. But since the homogeneous solutions are now bounded, this leaves the mathematical 
problem incompletely specified (i.e., ill posed) and the imposition of causality appears to be the most reasonable 
way of making the solution unique in this case. Our view is that the causality amounts to more than just imposing 
appropriate initial conditions and that it serves ( as demonstrated at the end of section 4) to insure the appropriate 
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cause- effect relation between the sound and its turbulent source when it is imposed on the Green’s function, as is 
done below. 

The present paper develops a causal solution to the relevant acoustic equations by using a vector Green’s 
function approach and assuming that the mean flow spread rate, say s, is small. This implies that all streamwise 
changes in that flow occur on the slow streamwise length scale ex u where X\ is the (suitably normalized) coordinate 
in the flow direction. The appropriate casual Green’s function consists of a component that decays to zero when the 
unsealed streamwise coordinate X\-X\ becomes large (where x{ corresponds to the source location) plus a 
component that becomes unbounded when the unsealed streamwise coordinate becomes large but decays to zero on 
the long (slow) streamwise length scale e (x\ - X\). The former component can be calculated by treating the slow 
variable SX\ as a parameter and using the locally parallel flow approximation to simplify the results. 

This approximation cannot, however, be used to determine the latter component (which corresponds to a 
superposition of linear instability waves; Briggs, 1964; Bers, 1975), since the resulting local solution would become 
unbounded and would not remain valid over the long streamwise length scale e (x\ — x/) on which this component 
evolves. The appropriate result is obtained by retaining 0 (e) terms in the Green’s function equations and using the 
method of multiple scales, the WKBJ method and matched asymptotic expansions to render the solution uniformly 
valid. 

Both components of the Green’s function act on the same source term and each is capable of producing acoustic 
radiation — even at subsonic Mach numbers. The first corresponds to the usual Lilley equation solution, but with 
slowly vaiying coefficients and with slightly modified source terms. The second is associated with linear instability 
waves, but is very different from conventional instability models since these waves are now continuously 
“generated” along the length of the jet and do not constitute separate sound sources. Their only role is to produce the 
appropriate Green’s function, and they may not even correspond to actual physical flow structures. Mani’s (1976) 
assertion that the instability waves generate the turbulence and should, therefore, be excluded is irrelevant here 
because their inclusion in the Green’s function only serves to produce the appropriate cause-effect relation between 
the sound and its turbulent source and, therefore, corresponds to a different aspect of their behavior. 

Each component of this function can be thought of as a filter that only allows certain elements of the source 
spectrum to reach the far field. The resulting acoustic spectrum is expected to exhibit a bi-model structure because 
the second Green’s function component only responds to frequencies in the range where the instabilities exhibit 
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spatial growth, which is somewhat lower than the frequency range selected by the first term. Preliminary 
calculations for a two dimensional mixing layer (Goldstein & Handler, 2003 ) suggest that the contribution of the 
second Green’s function component will be fairly small at subsonic Mach numbers, but that it should be quite large 
at supersonic speeds and dominate the first term at small angles to the downstream jet axis. Recent experimental 
results (e.g., see Tam, Golebiowski, and Seiner, 1996) show that the small angle spectra radiated by supersonic jets 
are quite different from those radiated at larger angles (say, at 90°) and even exhibit dissimilar frequency scalings 
(i.e., they scale with Helmholtz number as opposed to Strouhal number). Until now this curious behavior seems to 
have defied explanation, but our model calculations show that it can be attributed to the fact that each of these 
Green’s function components is the dominant carrier of the acoustic radiation at different angles to the jet axis. As 
noted above, the large differences between the present result and the usual non-causal Lilley’s equation solutions are 
directly related to the chaotic nature of the flow. 

The fundamental acoustic equations, which we refer to as the LNS equations, are introduced in section 2 and the 
Lilley’s equation as well as the equations for a steady but nonparallel base flow (i.e., the actual time- average flow) 
are obtained as special cases in section 3. A formal vector Green’s function solution to the general LNS equations is 
given in section 4 and a local causal, small- £ asymptotic approximation to this solution is constructed in section 5. 
The corresponding uniformly valid solution is obtained in section 6, and its far field expansion is worked out in 
section 7. This result is then used to derive an expression for the far field acoustic spectrum in section 8. The 
general formulas are applied to a round jet in section 9 and some simplifying approximations are introduced in 
section 10. Section 1 1 contains some qualitative comparisons between the numerical computations based on this 
simplified model and the available experimental data. 

2. The Fundamental Equations 

Goldstein (2002, hereafter referred to as I) showed that the Navier-Stokes equations, 


5 A 5 

— A H 

dt dxj 


= 0 , 


( 2 . 1 ) 


where the summation convention is being used, but with the Greek indices ranging from 1 to 5 (to represent the five 
first order equations), while the Latin indices i, j are restricted to the range 1,2,3, {A v } = { pv ; -,p/7 0 - p, p} , 

{ r v/} = {pv/v, + 5vP-<*ij, PVjho+q, -VjGy, pv, } , 
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(2.2) 


h = h + — v 2 
0 2 

denotes the stagnation enthalpy, h denotes the enthalpy, t denotes the time, A' = {xi,x 2 ,x 3 } are Cartesian coordinates, 
p denotes the pressure, p denotes the density, v = {vi,v 2 ,v 3 } is the fluid velocity, <jj is the viscous stress tensor, q, is 
the heat flux vector, { pv/,p/? 0 - p, p} is shorthand for { pvj,pv 2 ,pv 3 ,p/? 0 - p , pj etc., and the dependent variables are 
assumed to satisfy the ideal gas law: 


p=pRT, h = c p T , (2.3) 

with R = c p - c v being the gas constant, c p and c„ are the specific heats at constant pressure and volume, respectively, 
and T the absolute temperature, can be recast into the form of the linearized Navier-Stokes equations by dividing the 
dependent variables 


P = P + P'> P = P + P h = h + h', Vf =v i +v' i , (2.4) 

as well as the viscous stresses - // and heat flux q h into their ‘base flow’ components p,p,h,v i ,a i j, and q t , and into 
their ‘residual’ components p ’ ,p' ,b ! and q\, and requiring that the former satisfy the inhomogeneous 
Navier-Stokes equations 

^A v +^r v /=0, (2.5) 

at oxj 


along with an ideal gas law equation of state, 


h=c T = —— 
p Rp’ 


( 2 . 6 ) 
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where | A v | = | p v h ph 0 -p-H 0 ,p} ,{f vy - } = {pv/v,- +8 ij p-a ij -Ty, pVjh 0 +qj -v^y -Hj-vjHg, pry } 


h = /? + — v 2 
0 2 


(2.7) 


is the base flow stagnation enthalpy, and the ‘sources strengths’ 77, H a , and /? ; , which are assumed to be localized, 
can otherwise be specified arbitrarily. 

The residual variables are governed by the Linearized Navier-Stokes (LNS) equations 

V'v = ^ . (2.8) 

where 


{«v} ={p u'i,p'o, p'} 


(2.9) 


with 


P u'i = pv* 


(2.10) 


and 


P 


o = />'+(y-i)[^p v ' 2 



(2.11) 


is a five dimensional (non-linear) dependent variable vector, 

L V.V = + 5 h5 ) + K yiv > 


(2.12) 


with 


1 0X ■ 

^v-5 v ^---^5 v5+ (y-l) 
P St/ 


5v 


1 3f v 


K dx J 


J £ 1 vy 


J v4 — ~ 

P Sx, 


->4 




(2.13) 
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? <y =‘v ,-7 i/ _CT !/ ■ 


(2.14) 


D 0 denotes the linear operator 


d_ 
SX ; 


4 =— ,i = \i = 1,2,3 


(2.15) 


n d d . 

D = — + v • 

dt dxj J 


(2.16) 


and 0^ , , x ■ all equal to zero when p > 3, is the 5-dimensional linear Euler operator. The 5-dimensional source 


vector s M is given by 


5 , , j. / r 


where y = c p /c v is the specific heat ratio, the source strengths e' v are given by 

4 = -P vX - 4 + 5 /v ( Y - 1) I + H„ 

for v = 1,2.., 4 and zero otherwise and we have put 


+ < T'v> 


(2.17) 


(2.18) 




v ; =(y-l) /f + ^v ' 2 =c 2 '+^v 


4=(y-i)(4-^,) 


and 


(2.19) 

(2.20) 


4=-(y-! )(4-cr y/ v;) 


(2.21) 
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Equations (2.9) are easily converted into the usual convective form of the LNS equations by using the 5 th component 
of the base flow equation (2.5) to show that 


c " p/ - p a +v '&; f ' 

V J 


They are written out in the more familiar, but less compact, vector form in I. 


3. Lilley’s Equation and the Non-Parallel Mean Flow Result 

The base flow equations (2.5) reduce to the usual Euler equations when the arbitrary source strengths 
Tjj , H j, and H a , and the mean viscous stress ct (/ and heat flux tj l are set equal to zero. A general class of solutions 

to these equations, which quite conveniently provides a good approximation to the actual mean flow field in a jet, is 
the unidirectional transversely sheared mean flow 


Vj = 5 n U ( .r 2 , x 3 ) , p = constant, p = p (x 2 , x 3 ) . 


The 5 th component of the LNS equation then decouples from the remaining four components, which now 
become the inhomogeneous compressible Rayleigh equations (Betchov and Criminale 1967) 


- D 0 u'i , dU dp' d , 

P - 2 - L +8 n U: +-^ = 4 

Dt J dx , dx i dx , 9 

V J J 1 J 


D 0 p' 0 _du) de\j 8U 


where DJDt reduces to the usual convective derivative 


Dt dt dx x 


(3.4) 



It is now well known (see chapter 1 of Goldstein 1976), that the velocity-like variable a ■ can be eliminated between 


these equations (by taking the divergence of the first equation and the convective derivative of the second, 
subtracting the results and then using the first equation to eliminate the velocity fluctuation on the left-hand side) to 
obtain the inhomogeneous Pridmore-Brown (1958 ) equation 


, ' D » 
Lp ~ — 

0 Dt 


f 


de'j 


d Pz 
— c 

y 9x i dX ] J 


du 2 g 2 g 1 2 4 

dXj dxfiXj 


.K 

Dt 2 


de\ 


4 j 


dU 


dx. 




(3.5) 


where 


L^ 

Dt 


2 


d Ui d D ( 
dx: dx, Dt 2 


_2—— c 2 — 
dx, dx j dx. 


is the variable-density Pridmore-Brown (1958) operator, and 


c 2 =yp/p(x 2 ,x i ) 


(3.6) 


(3.7) 


is the square of the mean flow sound speed. 

This is a form of Lilley’s equation that was derived in I. As noted in the Introduction, it possesses homogeneous 
solutions corresponding to spatially growing instability waves on the base flow (3.1) (Betchov and Criminale, 1967). 
It also possesses homogeneous solutions corresponding to the continuous spectrum. The complete solution to this 
equation consists of a particular solution plus these homogeneous contributions, but the result is meaningless 
because the instability waves become unbounded (infinite) far downstream in the flow. The usual resolution to this 
dilemma is to completely neglect the instability wave contribution. Unfortunately, this causes the solution to be non- 
causal which, as noted in the introduction, can lead to erroneous results. 

A better resolution is obtained by choosing the base flow in (2.8) to be the actual mean flow of the jet. The over 
bars on the dependent base flow variables then denote the time average 

1 r 

•= lim —{•( x , t )dt (3.8) 

T^rao 2.1 J 
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where the dot is a place holder for p, v„ p, and h, and 


(p*/p) 


(3.9) 


denotes a Favre (1969) averaged quantity for all variables except h 0 , which is defined by (2.7). Notice that (2.6) is 
completely consistent with the overall ideal gas law (2.3) when the tilde is defined in this fashion. 

The time derivatives drop out of the base flow equations (2.5) and the source strengths are given by 

= -ptk (3.10) 


Ho = 2 T “- 


(3.11) 


The base flow equations are now the ordinary Reynolds-averaged Navier-Stokes (RANS) equations, which do 
not, of course, form a closed system. They are usually closed by assuming some sort of model relating the source 
terms to the mean flow variables Vj,p,p,and h and their derivatives, such as a Boussinesq model (Speziale 1991; 
Speziale and So 1998) for the Reynolds stresses and a similar model for H • . 

The most important difference between these results and the parallel flow result is that the homogeneous 
solutions to the LNS equations, which now correspond to instability waves growing and then decaying on the 
diverging non-parallel base flow, will always remain bounded. This eliminates the paradox alluded to above and the 
corresponding LNS equations can be used to calculate the radiated sound. The exact solution for a jet emanating 
from a nozzle consists of a particular solution that is driven by the sources (i.e., it satisfies causality and therefore 
the appropriate initial and far-field boundary conditions) and a homogeneous solution that is determined by the 
details of the nozzle geometry. The latter solution, which would be expected to vary considerably with experimental 
conditions and, therefore, be of only limited value as a predictive tool, is also somewhat inconsistent in that the 
instability waves must eventually become nonlinear in the actual flow while being forced to remain linear in the 
acoustic analogy framework. This solution has already been considered by Tam and Morris (1980), Tam and Burton 
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(1984) and others and only appears to be important at the high supersonic Mach numbers and/ or for jets at high 
temperature ratios (Mani, 1976; Tam 1995 ) where Mach wave radiation becomes dominant. We therefore restrict 
our attention to relatively low Mach number supersonic flows ( say, M < 1.5 ) and consider the former solution, 
which can be written in terms of the causal vector Green’s function for these equations. The general result for an 
arbitrary base flow is given in the next section. 

4. Formal Green’s Function Solution for the LNS Equations 

The particular solution can, as noted above, be expressed in terms of the vector Green’s function (Morse and 
Feshbach 1953, pp. 878-886) g m (vr,t|vr',t') , which satisfies 

L hv8vo =^ a 8(x-x')8(t-t'), (4.1) 

together with the causality condition 

g m (x,t\x',t') = 0 for t <t' , (4.2) 

and leads to the following formula for the dependent variable vector u v 


m v (x,t) = \ y J g V[l (x,t|x',t')^ ( x 'd')dx'dt ( , ( 4 -3) 

v 

where the symbol V denotes integration over all space. 

The derivatives acting on the source strengths el in (2.17) can be transferred to the Green’s function to obtain 


My ( x ,t) = “J ° J ~l vj / w dx'dt ' , 

V 


(4.4) 


where 
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(4.5) 


/ . \ d dv 

y vjtl ( x ,t\x\t’)^—g^-(y-\)— gvA . 

From the acoustics perspective, the primary interest is in the 4th (i.e., the pressure-like) component of (4.3), 
which, in view of (2. 14),(2. 1 8) to (2.21), (3.10) and(3.11) can be written as 


P'o = J m J Y jn {x,t\ *'/) t h/ (x',t')dx'dt' 

V 


(4.6) 


where 


d y - 1 ~ dg 4l , . 

' dx) 8 ^ 2 ^ dx\ + ^ Y ^ 


8x'j 


Y-l g 

2 ^ dx\ 


8 44 


(4.7) 


and 



(4.8) 


when the bulk viscosity is zero and the base flow is taken to be the actual mean flow in the jet. So in the inviscid 
limit, which is of primary interest here, Tj l1 is just a generalized four-dimensional fluctuating Reynolds stress and 

equation (4.6) therefore provides a direct linear relation between this quantity and the far -field pressure fluctuation 
(recall that p' a reduces to the latter in the far field). Notice that the causality condition (4.2), which requires that the 
Green’s function vanish when t is less than the emission time t' and does not depend on there being some initial time 
before all the fluctuations are “switched on”, merely serves to insure that all of the sound is generated by x • and 

that none of it comes from extraneous sources. 
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5. The Local Green’s Function Solution for a slowly Diverging Mean Flow 

The exact solution to (2.8) would require extensive numerical computation. But it is reasonable to seek an 
(approximate) asymptotic solution in the limit of small jet spread rate, say e, since most jet flows are nearly 
parallel — especially those of technological interest (see figure 1). In this section, we obtain an appropriate local 
approximation for this result. To this end, we suppose that all lengths have been normalized by some characteristic 
cross flow dimension of the jet (say its diameter, in the case of a round jet) and all velocities by some appropriate 
characteristic streamwise velocity with similar obvious normalization for the density, pressure and temperature. 

Then the mean flow velocities should expand like 

Vj = fy(A,jc ± ) + st/ (1) (X,jt ± ) + ... (5.1) 

v ± =sV(X,x L ) + s 2 V {l) (X,x ± ) + ... (5.2) 

where 

X = ex 1 (5.3) 

denotes the slow streamwise variable and 

*L={*2*3 }> ={^2^3} (5.4) 

denote cross flow variables. The expansions for the remaining mean flow variables are given in Appendix A, where 

it is shown that the mean flow is determined ( at lowest order ) by (A.. 7) to (A.. 12) in the inviscid limit 

where a (/ = q. = 0 -which is clearly appropriate here, because the mean Reynolds stresses are much larger than the 

viscous stresses in most cases. 

These results imply that the vector Green’s function g va and the linear operator L tlv in (4.1) must expand like 

g TO (x, x ,f|jc ; ', x ',/') = g^ ) +sgW +... (5.5) 

and 

V =4$ +eZ $+— ( 5 -6) 
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Substituting (5.5) and (5.6), along with the mean flow variable expansions, into (4.1) shows thatgj,^ is determined 
by 


+^-8 { S +8 a ^-g% ] = 8 ia S(x-x')S(t-t'), 


Dt 


dx, 


dxj ' 


(5.7a) 


= ^g ( 2 +^-(co 2 4" } ) = 5 4a 5 (*-*') 5 ( / -0» 


(5.7b) 






Dt 


dxj 


(5.7d) 


for ct = 1, . . .4, where Cq = , with P and R defined in (A. 1) and (A. 2) , and DJDt is now given by equation 

(3.4). The first order perturbations are given in appendix B. 

We begin by constructing a local solution that is valid in the vicinity of x 1 = x\ , but, as indicated in the 
introduction, does not remain valid on the long streamwise length scale X- X = X- sx{ and can therefore not be 

used in (4.6) to calculate the distant sound field. This solution is corrected in the next subsection to obtain a 
uniformly valid result that overcomes this difficulty. 

The left sides of equations (5.7a,b) are essentially the same as those of the parallel flow equations (3.2) and 
(3.3) and the velocity-like variables can again be eliminated to obtain the following single equation for g^ 4 J 


-2^ 0 2 ^8 ia 8(x-x')8(t-t')-(y-l)-^8 4a 8(x-x')8(t-t'), 

uXj uX\ Ut 

D2 O2 

where L is given by equation (3.6), but with c replaced by Cq . 
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Since the coefficients in these equations are independent of time and only depend on X\ through the slow 
variable X, a local ( causal ) solution can be obtained in the usual way by taking the double Fourier-Laplace transform 


(n) 

S\J = 


n 




m,k,X;x j 


f 

x dk dm for n = 0,1 

■ L , 


(5.9) 


where causality is enforced by requiring that the co integration be over the Bromowitz contour, i.e. the contour 
(shown in figure 2) that lies above all the singularities of g[°J . Appendix E shows that the 4 th component Fourrier- 
Laplace transform, g[°J ,can be expressed in terms of a single symmetric function 
G a (cb ,k,X;x ± |jc x ) = G 0 [m,k,X-,x' L |jc x ) that satisfies 


- _8(x L -x' L ) 


4 4 = 


(2n)- 


(5.10) 


where C k denotes the Rayleigh operator 


02 




4 = 


fai (kU-mf dxj ( kU-m )' 


- + 1 -- 


kC, 


7=2,3 


(5.11) 


It, therefore, follows from (E.3) that the 4th component of the lowest order Green’s function becomes 


(o) 

g\J = ■ 


n 


Co (4 


) 


\jiU (*^j_ ) — coj 


- d ; 


j[k(x, xf) to (< ( ')] G 0 ( | jci ) c/co 


(5.12) 


where 
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D' = 


, for ct = j = 1,2,3 

dx) 


-(y- 1 ) 

?>(*i) - 


m + U(x' L )— 


dx[ 


for ct = 4 


(5.13) 


Steady state solutions can only exist if the Laplace inversion contour can be continuously deformed onto the real 

00 

axis (otherwise J e ,k ^ x ~ x ^G^J dk would possess singularities in the upper half co-plane, which would cause g^j to 

-00 

grow without bound as t — > oo with x fixed). But jet flows are usually inviscidly unstable which means that G 0 must 

possess (usually simple) poles in the upper half /c-planc that cross the real /c-axis during this deformation. The k- 
integration contour c k must obviously be deformed to lie below these poles (as shown in figure 3) in order to obtain 
a continuous result in the co-plane (Briggs, 1964; Bers, 1975). 

The poles correspond to the infinitely many discrete eigenvalues, say k = k„ (to JO of the two-dimensional 
Rayleigh operator (5.11), but only a finite number of them cross into the lower half-plane (since there are only a 
finite number of unstable modes; Wundrow, 1996) when the cross section of the jet is finite. The integral over c k can 
then be decomposed into an integral over the real £-axis plus a finite number of contour integrals which can be 
evaluated as the sum of the residues of the poles that have crossed the real axis. 

Assuming simple poles, G a must behave like 


t T/t 

x ± x ± ,co,X 

k~ k„ 

where T„ is bounded at k = k„. Substituting this into (5.10), multiplying by an arbitrary function /(xj and integrating 
over all x± shows that 


as k — > k 


(5.14) 



J/(*±)4 r « dx l 


^ as ^ > k„ . 

( 2 ") 


(5.15) 
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But since /is arbitrary, this implies that £ K Y n = 0 and therefore that Y n = T n (a/,©, (jc ± ) , where (p„ is an 

eigenfunction corresponding to the eigenvalue k„. The corresponding spectrum may be degenerate with more than 
one eigenfunction corresponding to any given k„. It then follows from the symmetry property of G a that Y n must be 
of the form 


r„ = 




no sum on 


n ) 


(5.16) 


where A' n and ay = Uji are only determined when equation (5.10) is solved for G a . 
The causal Green’s function (5.12) can therefore be written as 


OO CO 2 / / \ 

JJ = - J J ( x± \ x '^dcodk 

^0 -00 [AT/ (a/ )-(£>]■ 


S* 


~ 2m L J r ~ ~ — - °a eL J( p m( x ±) dl 

» -o[k„C/(a:_l)-coJ A„ 


(5.17) 


where the sum is over the finite set of eigenvalues (which may be empty, i.e. the sum may be zero) that cross the 
real /-axis as the Laplace inversion contour approaches the real co-axis. 

6. The Uniformly Valid Green’s Function 

The solution (5. 17) consists of the sum of two terms, the first of which remains bounded, is a uniformly valid 
approximation to the true, non-parallel, result (or, at least, partially so, see end of section), and corresponds to the 
usual Lilley’s equation solutions that appear in the literature. But the second term grows without bound as X\ - x[ 
becomes large (since Im k„ < 0) and therefore becomes invalid on the long streamwise length scale X - X . In this 
section more or less standard perturbation methods are used in seriatim to render this term and, therefore, the entire 
Green’s function (5.17) uniformly valid everywhere in the flow. 

It can be made uniformly valid in the region x ± -x'j_ = 0(l) ,X-JC = 0(1) occupied by the jet by using the 
method of multiple scales ( see, for example, Crighton & Gaster, 1976) to obtain 
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( 6 . 1 ) 


2nd term in (5.17) — » -2ni^ j J 

n -oo 


S(*l)4,/W x ')%A x ±) a u D , 

_k„C/(^)-co] 2 a;(co,X') 


ip n j(x' ± )dm. 


where all quantities are evaluated at W, except where indicated and the amplitude function A nl satisfies the initial 
condition 


A nl (X |X')->1 asX^X' 


( 6 . 2 ) 


and is otherwise determined by imposing a solvability requirement on the first order solution which, as shown 
in appendix C, implies that 


!,( 0 ) 


c'A. 


dX 


nl + h\g A nl = 0, (no sum on n , /) , 


(6.3) 


where h. 


(o) 


and are given by equations (C.9) to (C.14). The solution to equation (6.3) that satisfies (6.2) is 


-\K,(x)dx 

A,i = e *’ ( 6 - 4 ) 


where 



(6.5) 


While (6.1) provides a uniformly valid approximation to the second g^j term within the jet, the result is still 

not valid at large transverse distances because the solution tpj^ of (C.3), which is a generalization of the second- 
order mode shape equation in the linear instability wave problem considered by Tam and Burton (1984b), decays 
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more slowly with |jc ± as |x ± | — > oo than the first order solution tp 4jl/ (see also Tam and Morris, 1980). Tam and 

Burton consider only a round jet, but the large |jc ± | asymptotic behavior of the present, more general, result is the 

same and their demonstration applies to the current situation as well. They also obtain an appropriate outer solution 
that remains valid at large transverse distances by using a Fourier transform approach but in the present paper, we 
adopt a somewhat different approach, based on the WKBJ solution, to the reduced wave equation in order to deal 
with the more general (non-axisymmetric) flow being considered herein. 

To this end we need only consider the factor 


. x 

— f YL n (to,X)dX 

J ! = A nl (x\X')e £jr (fcrfK) ( 6 - 6 ) 

in (6.1). Once its far field expansion is known, the result can be substituted back into this formula to obtain the 
appropriate far field behavior of the corresponding component of the Green’s function (5.17). Since 


/ \ 

2 


f \ 

2 1 

[ £® 


vi+| 

CO 

(N 

1 




v c oo y 



as r ± = jcJ — > oo , 


(6.7) 


where 


2 i s 5 i r 

V l=— T-r ± — + — - T , 
r L or L dr L r[ 5cp 


( 6 . 8 ) 


(p is the azimuthal coordinate and c x is the free stream sound speed, it follows that <p„/ must behave like 


1 , . \ Jk 1 --(oy'e. ) r 


as r ± — > oo , 


(6.9) 


which means that 
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1 


as r ± — > oo , 


( 6 . 10 ) 


J, 


A(X |X')<t((p)e 


1 X 

-\*(X)dX+ 


where we have, for simplicity, dropped the subscripts n and /. 

Since qv arises from the continuous deformation and uniform extension of the original Green’s function 
formula (5.9), the branch cuts for the square root in (6.10) must be consistent with that result (see equation (7.2) 

below), which means that -(ro/c^)" must be positive when k 2 > (ro/c*) 2 and negative imaginary when k 2 < 

(cb/Coo) 2 . This only occurs when the portion of the branch cut emanating from k = © /c® extends to infinity in the 1st 
quadrant and the portion emanating from k = -co/c.,. extends to infinity in the 3rd quadrant. Their detailed shapes are 
arbitrary at this stage of the analysis. 

The result may not remain bounded as r± — > oo when Im k (X) > 0 (i.e., when the instability wave is damped, see 
Tam and Burton, 1984) but, as shown in section 7, it still produces an outgoing wave in the far field with this choice 
of branch. The only requirements are that the inner and outer solutions match and that the boundary condition at 
oo be satisfied. The outer solution eventually converts the unbounded inner solution into a propagating or damped 
disturbance. 

The outer solution, say J a , that is valid in the region 

X 2 +Rl= 0(1), (6.11) 


where 


R_ i = K'l , 


( 6 . 12 ) 


and matches on to this result when R ± — » 0 with X> 0 must satisfy the reduced wave equation 
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v d 2 + 


1 8 _ 
R 2 l Sep 2 


2 > 


\J n + 


f A 2 

05 


J 0 =<>, 


(6.13) 


where 


v d 2 = 


—R, — — 

S7? ± 5i? ± 


dX z 


(6.14) 


and must be of WKBJ form 


J n =^^ 0(XA) o((p)4(x,i? 1 )+... 


v/^i 


(6.15) 


Substituting this into (6.13) and equating like powers ofg shows that 


0^+© 2 ± -(cq/ C oo ) 2 =O 


(6.16) 


and 


0, 


S 

8R, 


+ ©, 


dX 




7 0 z © = 0 , 


(6.17) 


where the subscripts denote partial derivatives with respect to the indicated variables. 
The solution to the Eikenal equation (6.16) that matches with (6.10) is 


0= J k (Z)dZ + - 


X' 


(to /gpo ) 2 R 1 
-(“/Coo f 


IlK' 


(6.18) 
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where 




(6.19) 


=X-a R ± 


( 6 . 20 ) 


and 



( 6 . 21 ) 


which can be verified by differentiating (6.18) to (6.21) to obtain 


©v =k 


( 6 . 22 ) 


and 


0 


R ± 



(6.23) 


Expanding (6.18) in a Taylor series shows that 


0= j K(X)dX + i^K 2 (X)-((o/c x fR ± +0{ y Rl^j asR ± ^0, 

X' 


(6.24) 


which, in turn, shows that (6. 18) matches with the phase factor in (6.10) as R± — > 0. 

It is worth noting that a satisfies the (in general complex) inviscid Burger’s equation (Burgers, 1948; Pierce, 
1981, pp. 588-591) 
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subject to the (initial) boundary condition 


a(X,0) = a(X) for X > 0 


(6.26) 


The branch of the square root in (6.18) must agree with the choice of branch in (6.9) and (6.10) and must otherwise 
represent a continuous function of X and R to the maximum possible extent. The Burger’s equation solution, which 
always develops a singularity when its dependent variable is real (Pierce, 1981, pp. 588-591), may also develop a 
singularity for complex a when the mean flow Mach number is sufficiently large — corresponding to the sonic point 
singularity, k = ro/c^ , in (6.19). The resulting solution would then be discontinuous across a line connecting this 
point to infinity, which can be made to coincide with the branch cut for the square root in (6.21) and thereby 
minimize the discontinuities in a . These discontinuities correspond to the caustics that arise in the geometric (i.e., 
high frequency) acoustics (Avila and Keller, 1963; Pierce, 1981, pp. 381-384) and signify a local breakdown in the 
WKBJ solution (6.15), which can be eliminated by constructing a local inner solution to the reduced wave equation 
(6.13), but we do not pursue this issue further. 

The amplitude equation (6.17) can be written more simply as 


d 

dR ± 


d® d 2 d® 

+ — A: — 

dR L dX 0 dX 


(6.27) 


Inserting (6.22) and (6.23) and using (6.21) to rearrange the result shows that 


— ad 2 +-^—Al = 0 . 

8X 0 dR ± 


(6.28) 


But differentiating (6.21) with respect to X shows that 



(6.29) 


5a _ a'(^ 0 ) 

8X ~ \ + a\X 0 )R L ’ 


where the prime denotes differentiation with respect to the entire argument X 0 . It is now easy to verify that the 
solution to (6.28) (and therefore (6.27)) that matches with the amplitude of (6.10) as R± — > 0 is 


a 2 (x,r ± ) 


A 2 {Xq\X') 

[l + a'(X 0 )7Jj_] ' 


(6.30) 


A uniformly valid composite expansion, say J Io for / that reduces to (6.6) in the inner region and (6.15) in the outer 
region (6. 11) is (Van Dyke, 1975) 


A n A X o\ X ') 


-®„(*A) 

e z 


<P„/(*l) 


(6.31) 


where 


©„=©„- i ^ 2 ( x )-Hc x ) 2 R 


(6.32) 


Substituting this into (6.1), and using the result to replace the second term in (5.17), yields the following expression 
for the Green’s function 


gia = j e ^(?4^(a:|a:';co)/co (6.33) 


where 
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CO D2 / t \ 

' L -,(D,k,X,X') = -\ C ° [Xl) D 'J k ^G 0 { Xl _\x' L )dk 


, \jdJ (x^ ) - 00] 


0 V"' C o (•*'±)^n/ {Xg \x )fl/y(p«/ (j£j_) i@(X,R)/e 1 , \ 

2 ni L~ r : — : — ^2 7 = D o e V( p„/K) 

» [k„c/(jc!)-co] a;^/i+< 


(6.34) 


the sum is over the finite number (which may be zero) of unstable eigenvalues at A" and frequency co, G 0 , A', , and 
a/j are determined by solving (5.10), 9,,/ and k„ are obtained by solving the eigenvalue problem for the Rayleigh 
operator (5.11) and ©„ ,X m and a are given by (6.18) through (6.21), and (6.32). 

Unfortunately, (6.34) is still not uniformly valid in X for all co .This is because the pole at k„ of the nth 
instability mode coalesces with its complex conjugate k* on the real £-axis when the slow stream-wise coordinate 
X of the source point .v' is equal to the neutral stability point X Nn , which causes the integrand in the first term of this 
result to become infinite there. The integral itself remains finite — but the local contribution from k* can be of the 
same order as the second term in (6.34) and must therefore be treated more carefully. To this end, we first separate 
out this contribution by noting that the original integral is the same as the integral over a contour that passes just 
above k* plus its residue 


□2 ( f 

M/Pj 

t ' 

*1,K„ 

t } 

) D > 

nj ( 

t * \ 


_<u0 


2 

a;,i 

(< 

) 


(6.35) 


at this point and that this latter quantity produces the dominant portion of the localized contribution. 

The numerical results show that A', (k„) and A'„ |k* j remain finite at this point and that A', |k* j = -A', (k„ ) 

there. The term (6.35)-which is itself a homogeneous solution to the governing equations-will therefore cancel the 
/7th term of the second member of (5.17) at the neutral frequency point when x l is also at this point. But this result is 
non-uniformly valid in X and does not cancel the appropriate member of the second term of (6.34) when X ^ X Nn 
even when r± = 0(1) (i.e., within the jet). It is however, possible to construct a uniformly valid result that completely 
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cancels the appropriate member of the second term in (6.34) when X Nn ( X' = X Nn ) by first constructing a 

homogeneous solution to the non-parallel Rayleigh’s equation that is uniformly valid within the jet (i.e., when r± = 
0(1)) and coincides with the local homogeneous solution (6.35) when X is near any of the neutral points X Nn . 

This can be done by replacing the factor (p n/ |x ± ,k x ^ in (6.35) by 

i 

A n\ x \X\ k* ) (p„ ; , < ) D’ a e L 

The result will be exponentially small compared to the second term in (6.34) when r ± = 0(1) andX < X Nn and will 
cancel the appropriate member of that term when X =X Nn . It can easily be extended to the exterior of the jet by 

f} K ,(w,x)dx e i®„(x m ,R,)/s 

replacing e * m by . 

V 1 + a n {X a )R L 

The resulting uniformly valid homogeneous solution to the non-parallel Rayleigh’s equation can be added to the 
uniformly valid result (5.34) to obtain a solution that remains valid even when.Y' =X Nm i.e. for all values of ro .It is 
important to note that the difference between this result and the original solution (6.34) will be exponentially small 
(i.e., asymptotically negligible) unless X is close to the neutral point X N „. 

7. Far Field Green’s Function 

The primary interest is in calculating the jet noise at large distances from the flow where the far field region 
where 


J Kl((d,X)dX+ j k j; (c o,X)dX 


sr = R — > oo 


(7.1) 


and the uniformly valid reduced Green’s function (6.34) takes on a much simpler form. The appropriate result is 
derived in this section. To this end, we note that, since G a behaves like 


-Jk 2 -(<o/c„y 


-Q 0 (cp,A:,co|jc^) as r ± — » oo , 


(7.2) 
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where the choice of branch has already been discussed in conjunction with (6.9) and (6.10), the integral in the first 
term of (6.3) can be readily evaluated by the method of stationary phase (Carrier, Crook, and Pearson, 1966, p. 274) 
with the stationary phase point occurring at 


cos Q-k 


smG/ijk 2 -(co/c x ) 2 


= 0 


(7.3) 


or 


k = 


03 Q 

— cost) 


(7.4) 


where 0 is the polar angle 


- r, R, 
tan 0 = — = — 
x l X 


(7.5) 


measured relative to the downstream jet axis, and 


R 2 =X 2 +R 2 i 


(7.6) 


is the radial coordinate. 

The far-field expansion of the second term of (6.34) is evaluated in Appendix D. The results, along with (6.9) and 
those of this section can now be inserted into (6.34) to show that 


? ( 0 ) 




-» — 


ff> 


sin0 


G Ha \x\x 


(%') + £el?Pk) 


(7.7) 
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where 


yjlni 


U2 / r \ 


co 


J B ( j 


[M(x^)cos9-l] 


J 7 -^ t -/COX.' COS0/c 


Go 


cp, — cos9,cd|x^ 
c„ 


(7.8) 


An) - 2ni °° KK/ K,» I X ') a IJ°nl (<p) n , J 0 :( e,X’)/s„ ( ^ 

“Ter ~ /- r . . -|2 I 72Z — ‘P/y 

VS [k/Z^D-Co] Kj<( a on) 


(7.9) 


x = x(9,cp) denotes the unit vector in the x-direction. 


a 0 

©« ( 0 ,-T)= J k„ (Z)dZ- — a o; , cos 9 

V' C '00 


(7.10) 


a 0 „ is determined implicitly by (D.2), (D.3), and (D.7) once the local growth rate K n (X,co ) is determined from the 
eigenvalue problem (C.8), and 


M(x' ± ) = U(x' ± )/c r 


(7.11) 


is the acoustic Mach number. Appendix D also shows that (6.15) exhibits the desired outgoing wave behavior at 
infinity. 


8. The Acoustic Spectra 

Jet noise theory is primarily concerned with predicting acoustic spectra, which correspond to the Fourier 
transform of mean square pressure correlation 
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(8.1) 


P 


2 


2 T 


T 


J P'o{ X ’ t )p'o( X ’ 


-T 


t + t Q )dt 


in the far field where r — » oo and p' Q , which is defined by (2.1 1), reduces to the pressure fluctuation. Here, T denotes 
a large time interval. A relatively simple formula for this entity will be obtained in this section. 

Inserting (4.6) into (8.1) and changing the integration variables to t\ = t — t" and x = t" — f, yields 


A- 2 = J JJ W^m{ x \ x 'd + t 0 -t')^ j]i (x\x",t-t")x ai {x',t')x VL j(x",t")dx' dx"dt' dt"dt 


-T -00 V 


(8.2) 


= }J JJy/a(- v |- v '^i+l 0 + T )Yyi- v k"^i) ^ iw ( x ',x"-x',x)dx'dx"dt l dx , 


-00 V 


where 


T cw( A: ' ;/ ‘ Jt '- T ) E ^ J T oi( X '’ t ') X w( X "’ t ' +X ) dt ' ( 8 - 3 ) 

-T 

is the density weighted 4th order two point time delayed turbulent velocity/total enthalpy correlation. Inserting (4.8), 
with a^y = 0 ,into this result shows that 


[x ;x -x 


> T )-<i(*'; 0 > 0 Ky(*'';°> 0 ) 


(8.4) 


where 


(x'-,x"-x',x) = ^;^ p(x',t')p(x",t' + x)v' a (x',t')v' i (x',t')v' il (x",t' + x)v'j(x",t' + x)dt' (8.5) 


-T 


and 


29 



( 8 . 6 ) 




are the 4th and 2nd order correlations of the fluctuating (i.e., turbulence) velocities or total enthalpy. Notice that 


f a/ 0,0) . (8.7) 


Equation (8.2) can be written more compactly as 


CO 

P 2 (x,t 0 ) = J JJy /a , ; (x|^,ri,t 0 + T) ^ {\’-,%x)dx’d\\dx , (8.8) 

-00 V 


where 


00 

Y/a/n = J y ;a (jc|jc',h + t 0 + x) yj*|*' + q,h) dt x , 

-oo 


(8.9) 


which can be expressed in terms of the Green’s function, accounts for the acoustic propagation and mean flow 
interaction effects. These results provide a direct linear relation between the mean square acoustic pressure and the 
4th order correlation of the turbulent fluctuations in the jet. Unfortunately, determining this latter quantity, either 
experimentally or numerically, is a very difficult task and it is, therefore, highly desirable to make the acoustic 
predictions as insensitive to its details as possible. In fact, this is one of the main reasons for using the moving media 
acoustic analogies rather than the simpler Lighthill formulation. This desensitization would certainly be enhanced if 
the variation of x . ; with q was much more rapid than that of the propagation factor y • ; , because the latter 

could then be treated as a constant relative to the q -integration and the result would then depend only on the 
integral of x • z over the separation variable q , which, in turn, depends only on the decay time x at any given 
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source point x' . This may be the case at lower Mach numbers, but the rjj -decay rate is expected to be much slower 
at the higher Mach numbers of technological interest. However, Lighthill (1952, 1954) pointed out that the 
streamwise decay rate should be much more rapid in a reference frame moving with the convection velocity, U„ of 
the turbulence and Ffowcs Williams (1963) showed that this idea can best be implemented by introducing a moving 
frame correlation, say 


*a jt xl = T q/W (*'’ t l’ T ) ( 8 - 10 ) 

where 

% = x\-i U c x , (8.11) 

before integrating the turbulence correlation over the separation vector. Introducing this into (8.8) and changing the 
integration variable to if (not to be confused with the c introduced in Appendix D) yields 


p 2 ( x do) = } jj Yjain( x \ x '& + i U c x,t 0 +x} x^,(x'-,^x)dx'd^d x. (8.12) 

-CO V 

This is a generalization of the result obtained by Ffowcs Williams (1963) in the Lighthill context. It turns out to be 
much simpler in frequency space when the observation point r is in the far field. 

To obtain the appropriate result, we introduce the far field acoustic spectrum 



-CO 


p 2 (x,t 0 )dt 0 . 


(8.13) 
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00 

note that the Fourier transform of a correlation, say J f(t + x)g( t) dx is just 2 jlF(o>)G*((o) (where capital letters 

-00 


denote Fourier transforms of the corresponding lower case letters, and the asterisk denotes complex conjugates), and 
take the Fourier transform to (8.12) to obtain 


Ao H*') = 271 J J r va(- v |- v ';®) r ^(- v |- v ' + ^ + l ' U c x\ay ™ Xa jn i{x'&,T:)d£, dx, (8.14) 

-ooV 


where the left hand side is defined by 



(8.15) 


and therefore denotes the acoustic spectrum at .v due to a unit volume of turbulence at x'. 
Now it is clear from (4.7), (6.33), and (7.7) to (7.9) that 


iar/c [— 

e ' ro • Q 
F = — — sin0 

^ r Vc. 


r j/a(*|*') + E r SS(*|*') 


(8.16) 


where 


x[ COS 0 

r Bj a =rBj a (x |*i) e C " =D JoVL G BYL (k\x'), 


(8.17) 


p( H ) _ 

1 Ijo - 




>:(0x) 


=^/c V G £ ) (*!*') 


(8.18) 


and (for the small spread rate being considered here) 
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provided 0 and M are sufficiently large. The result that is uniformly valid when the source point is in the vicinity of 


a neutral stability point is obtained by replacing 


_K„£/(*i)-«>] a ;(k„) 


D o e v " <P»y(*±»K n ) by 


[k„c/(jc;)-co] a;(k„) 


- Die 




_k*c/(jcI)-co] a;( k ;) 


-Die 


j <(co ,X)dX-KS>”(e,X Nl1 ) 




in (7.9). 

Since 0“ varies on the slow streamwise length scaled, which is much larger than the correlation length of the 
turbulence, we can account for its variation over the range of integration in (8.14) by expanding it in a Taylor series 
and using(6.22) to show that 


© 00 (e,X , + e(^ 1 +C/ c T)) = © 00 (e,X , )-eK(X')(^ 1 +C/ c T) + ... (8.20) 

with more than enough accuracy to evaluate (8.14). It therefore follows from (8.17) and (8.18) that 

r Bja (x\x' + $ + i U c x) = r Bja (x \x' + ^)e- ie>M ^ (8.21) 

rg(i|*' + l ; + / U c x) = T^{k\x' + %)e- iK - {X ' )U ^ (8.22) 


where M c = U c /c w denotes the convective Mach number. 

For simplicity, we only consider the case where a single mode is excited at a:' and ro. Other modes will, of 


course, be excited at lower frequencies and will have to be accounted for in order to predict the true acoustic 
spectrum. Then inserting these into (8.14) via (8.16) shows that 



(8.23) 


271 1 fO (* r nr 

/«,(*!*')= —) — sin0 J[ r ^(^|^O +r ^( i k')]L r W(- i 'l x,+ ^) // ^/(- v,; ^®( 1_M c cos0 ))' 


+^)Km (*'; ^ co - u c k* (x', ©)) 


for r » 1 , 


where 





(8.24) 


denotes the two point fourth- order turbulence spectrum relative to the moving frame and the second terms in the 
square brackets (i.e., the F //f7 terms ), which are identically zero when co is greater than the neutral frequency 

(where Im k(jc',<d) = 0), can be modified in the obvious way to obtain an equation that is uniformly valid at the 
neutral points X = X Nn . 

This result is, aside from the easily relaxed constraint of single mode excitation, nearly exact. The only 
significant approximation is that the mean flow spread, s, is small. But // 0/o/ is expected to vary much more rapidly 
with £ than Y ^ and (8.23) can therefore be approximated by 



— sin0[r^ a (*|*')+r j j ta (*|*')][r^(*|*')©* ;M/ (*';a)(i-iif c cos0)) 

C 00 


(8.25) 


where 


d’q/y {x';<a) = \ H ajv d (x';$,<o )<% (8.26) 

v 

is the single point 4th order turbulence spectrum at x' and the second terms in the square brackets must again be 
modified to make the result uniformly valid at the neutral point X = X Nn . 
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This approximation is also invalid when M c cos9 = 1 because the first spectral function in (8.25) does not go to 
zero as (fl — > oo in this case, which causes the integral of (8.25) with respect to ro and therefore the mean square 
pressure to become infinite. Ffowcs Williams (1963) showed that this difficulty could be overcome within the 

Lighthill (1952) context by replacing the Doppler factor (1 —M c cos9) by ^(l - M c cos 9)“ + (a M c )" where a is a 

relatively small constant. This result does not remain valid in the present more general context, but our interest here 
is in the case where M c < 1 

9. A Specific Example 

To make these results more concrete, we consider the case of an axisymmetric mean flow, i.e., a round jet, 
where 


U=U(r L ,X), 


%=%(r ± ,X) 


(9.1) 


and V is purely radial with magnitude V(r ± Jf). Then the mean flow equations (A. 7), (A. 8), and (A. 10) decouple 
from (A.9), i.e. they reduce to the usual boundary -layer equations, but with DJDt now given by 


Al 

Dt 


dX r L dr ± 


r L V 


(9.2) 


in place of (A. 1 1) and S given by 



in place of (A. 12). The radial pressure equation (A.6) reduces to 



(9.3) 


(9.4) 
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which can be used to eliminate P in (9.3), and the Reynolds stress components are defined in the obvious way. 


However, the - pj term 


is fairly small and is usually considered to be negligible (Pope, 2000, p. 114). 


This does not cause any in consistency in the asymptotic results and we shall suppose this approximation to be valid. 
The Rayleigh operator (5.11) now becomes 


1 d r^cl d 

r L dr± (kU-af dr ± 


R 2 


(kU - co) 2 


k 2 


1 d L 

r 2 5(p 2 


+ 1 


(9.5) 


and G a must therefore be of the form 


= I 


o m(cf>-<p') „(h) 


(fi | r i) > 


(9.6) 


where 


£ nk G {n) = 



(9.7) 


and 


It follows that 


1 d r ± cl d 
r ± dr L ( kU-a )f dr L 


1 -- 


c 2 


(, kU-(o ) 




(9.8) 


(«) = w l X) W, k Wn ] 

K{k) 


G K ’ = 


no sum on 


n ) 


(9.9) 
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for r ± > r\ ( and a similar equation with r ± and r[ interchanged when r ± < r ' ), where 


£ n k w l!\r L ,k) = 0 for/ = 1,2 


(9.10) 


./O 



as r, — > oo 


(9.11) 


,/ 2 ) 


— > r 


as r, — > 0 


(9.12) 


and 


A„ S -(2ji) 2 wJ { kU - u >) 2 , 


(9.13) 


where 


W „ 


(l) d ( 2 ) ( 2 ) d (i) 


dr , 


dr , 


(9.14) 


is the Wronskian of w! 


(i) 


and w! 


,( 2 ) 


, is independent of r ± . 


Since k„ depends on n 2 , the eigenvalue problem for C k is doubly degenerate for each n 1 (corresponding to ± rt) 


9,,/ -> <p,± = 


,( 2 ) 


( r l| K n)< 


+in ( p 


(9.15) 


and a/j is equal to 1 for I # j and zero otherwise. The constant A', that appears in (5.17), (6.1), (6.34), and (7.9) is 


now the derivative of (9.13) with respect to k, 


A' 


lim 

k—>K n 


d K[ k ) 

dk 


(9.16) 
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evaluated at k„ and it follows from (7.2), (9.6), (9.9), and (9.1 1) that <D, V and the function Q a in (7.2) are given by 


O . = e m>9 

nj 


(9.17) 


and 


So = X 


oo /«(cp— (p') 

e ,( 2 ) 




(9.18) 


All other quantities in (7.8) and (7.9) can be calculated from the preceding equations once the eigenvalue and 
boundary value problems for the operators ^nk have been solved. 

10. Some Approximate Results 

Equation (8.25) is very general and relatively exact but requires a great deal of modeling and/or computation is 
required to evaluate all of its terms. It is therefore worthwhile to introduce some approximations that can be used to 
gain insight into the importance of the various physical effects described by this result, which will be done in this 
section. 

Considerable simplification can be achieved by assuming that 




(10.1) 


which implies, among other things, that h' 0 1 0 . It follows from (7.8), (7.9), (8.17), (8.18), and (8.19) that 
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(10.2) 
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and, again assuming only a single excited mode, 
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(10.3) 


when o) is less than the neutral frequency at X' ( Im k(<j),X ) = 0 ) and equal to zero otherwise. Then, using (10.1) and 
neglecting the cross coupling terms in (8.25) yields 
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where the factor 
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to make this result uniformly valid in the vicinity of the neutral points 7 = X N . 

11. Some Numerical Results 

In this section, we carry out some numerical computations based on (10.4). Low Reynolds number direct 
numerical simulations (Freund, 2002) suggest that <D 0 can be fairly well represented by a Gaussian distribution, say 


O,, = ■ 


/3 

l c a c 


CD P 


f ^ 2\ 

PA 

4n 


— (co/2o: 


(ii.i) 


where p s is a constant representing a characteristic density, to v is a frequency (inverse time) scale of the turbulence, 
l s oc X' is a characteristic length scale, u s oc j/^., is a characteristic turbulent velocity, and a s is an empirically 
determined scale function. 

The experiments show that the mean velocity is reasonably well represented by (Tam and Burton, 1984) 

U/U cl = U(j\)/U cl (X') , (11.2) 


with 


n = r ± /X' . 


(11.3) 


over much of the jet. 

Pope (2000, p. 1 16) shows that this type of solution is consistent with the mean flow equations (A. 7) - (A. 10), 
for the axisymmetric case being considered here. Solutions to the Rayleigh stability equation for this type of profile 
depend on the frequency and the local streamwise coordinate only as the product a> = to A ' . 

We expect U c j to be nearly constant up to the end of the potential core (see figs. 5 and 6 of Tam and Burton, 
1984) and take U to be of the form 
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c/ ( 1 l) = ^[ 1_tanh ( 1 l“ 1 l c )] 


(11.4) 


in this region. q r = constant is the centerline o f the shear layer and the hyperbolic tangent changes rapidly enough 

to ensure that U (0) is reasonably close to unity. The mean sound speed profile is taken to be uniform and equal to 

its ambient value. While the present formulation is intended to minimize the dependence of the results on the details 
of the turbulence, they are expected to remain strongly source dependent and the numerical computations based on 
this simple flow model should only provide a qualitative indication of the expected experimental results. 

The imaginary part of the phase factor 0 00 in the second term in (10.4) arises from the WKBJ solution and 
accounts for the initial transverse decay of the instability before it turns into a propagating disturbance. This term 
will be negligibly small (compared to the first) when Im ©°° is 0(1) and greater than zero, and very large (i.e., it will 
be dominant) when Im ©°° is 0(1) and less than zero, since the spread rate e is very small (of the order of 0.1). 

F igure 4 is a plot of Im ©“ as a function of 0 for various values of the Mach number M and source frequency a> for 
the first three azimuthal instability modes. It shows that Im ©°° is always positive at 0 = 90° but can cross the axis 
and become negative when 0 is sufficiently small and M is sufficiently large. The results suggest that the 
n=l azimuthal mode will make the largest contribution to the instability wave component of the Green’s function, 
and that its magnitude will be greater at lower frequencies. 

The curves are truncated at certain relatively small values of 0, because a does not exhibit the asymptotic 
behavior (D.l) and (D.2) beyond these points and the second term becomes non-radiating ( i.e., evanescent) there. 
The 90° spectral shape is therefore produced by the first term in (10.4), which as noted above, corresponds to the 
usual Lilley equation solutions that appear in the literature and is primarily determined by the factor 

\ 4 T 03 V 

— ® 0 (a)(l-M c cos0)) = — 0 0 ( oo) at that angle. It has a similar shape at all other angles but its 

'-'co y V '" 00 y 

magnitude increases and its frequency scale is stretched when 0 < 90° — implying that the spectral peak moves to 
higher frequencies as 0 — » 0, which is not consistent with experimental observations (e.g., Tam et al., 1996; also see 
fig. 7(a) of Freund, 2002). 
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Some typical measurements are shown in figure 5 taken from Noram (1994). The jet Mach number is 
approximately 1.5 and the jet is fully expanded- which means that there is very little shock-associated noise in this 
case. Notice that the peak of the small angle spectrum is at a much lower frequency than that at 90°. When Lilley’s 
equation was first introduced (Lilley, 1 974) the expectation was that this behavior would be produced by the 
refraction effects that are included in the conventional bounded solution (which corresponds to the first term in 
(10.4))- the argument being that they would cut off the high-frequency part of the spectrum at small angles to the 
downstream axis and thereby cause an apparent shift to lower frequencies. But our numerical results (as well as 
those of Khavaran and Bridges, 2004) show that all the relevant frequencies are equally reduced by refraction at the 
higher Mach numbers being considered here and the net result is an increase in the peak frequency of this 
spectrum — which is just the opposite of what is shown in figure 5. But figure 4 shows that the second term can 
become dominant at these angles when M is sufficiently large. Its spectral shape is relatively independent of angle 
and has a much narrower width than <D 0 (a>) due to the relatively narrowband of unstable frequencies at X . 

The argument, a>-[/ f K(Jf',a>) , of <D 0 tends to be close to zero over most of this range (where the second term 

is dominant) — which means that the resulting acoustic spectrum only depends on the low-frequency component of 
the turbulence and is therefore relatively independent of its spectral shape. This is in contrast to the spectral shape 
produced by the first term which is primarily determined by the turbulence spectrum <D 0 . The second term is 
expected to underpredict the low frequency component of the spectrum where the appropriate Green’s function 
would involve more than a single instability wave mode. 

Figure 6 is a plot of the spectral shapes computed from (10.4) forM= 1.5 with individual contributions from 
the n=0 and n= 1 azimuthal instability modes (but neglecting the cross coupling terms) at the polar angles 
corresponding to the experimental data in figure 5. We have taken M c = 0.65 M, fj c = 4 ,and a> s = 0.3 in all of these 

calculations. The spread rate is calculated from s = (0.165 - 0.045 M 2 )/ -Jn (Lau, Morris, and Fisher, 1976). 

The computed results, which are for a ring source at q' = 2.5 , are seen to exhibit the general characteristics of 
the data shown in figure 4: namely the spectral peak at 0 = 30° is about a factor of ten larger, and occurs at about 
half the frequency of that at 0 = 90°. Each of the two terms in (10.4) produce a distinct lobe in the small angle 
spectrum, with the low frequency lobe resulting from the second of these. The abrupt change in slope can probably 
be attributed to our neglect of the interaction terms. The long-dashed curve is our estimate of the effect of retaining 
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those terms in the spectral calculations. Of course, the actual acoustic spectrum produced by the jet is a sum over 
these individual point contributions and will probably not exhibit distinct lobes — even when the interaction terms 
are neglected. The summation should also produce a smoother spectrum. 

The M= 1.5 spectral shapes shown in figure 6 are plotted at a number of additional polar angles in figure7. The 
instability wave contribution (i.e., the second term contribution) to the acoustic spectrum maintains the same shape 
but increases in magnitude as the polar angle decreases. This contribution disappears rather abruptly when 
8 becomes larger than about 50°. Figure 8 shows the spectral shapes for M= 1.3. 

A major paradox that seems to have eluded explanation is that, while the 90° spectral peaks scale very well with 
Strouhal number (Gaeta & Ahuja 2003) as predicted by Lighthill's theory, the small angle spectral peaks observed 
in high subsonic and low supersonic Mach number experiments scale much better with Flelmholtz number (Ahuja, 
1973 ; Tam et. ah, 1996) . Figure 9 shows that the 90° spectral peak is unchanged while the small angle spectral peak 
shifts to lower frequencies as the Mach number increases, which means that the 90° spectrum scales with Strouhal 
number but that the small angle spectrum does not. The latter is replotted against the Flelmholtz number in figure 10, 
which shows that the predicted scaling of the spectral peaks is now consistent with experiment. 

12. Concluding Remarks 

The Navier-Stokes equations, rewritten in the form of the linearized Navier-Stokes (LNS) equations with 
externally applied stress and energy flux sources, which we argue are the natural generalization of the Lighthill 
(1952) zero mean flow acoustic analogy, were solved by using a vector Green’s function approach and assuming 
that the spread rate e of the mean flow about which the equations are linearized is a small parameter. The relevant 
solution satisfies causality and the resulting Green’s function consists of two components one of which involves 
spatially growing instability waves. The solution is therefore very different from the parallel flow results even 
though e is very small. Numerical results were obtained for a simplified model of a round jet. They show that the 
small angle spectrum is narrower and has a lower peak frequency than the 90° spectrum which is consistent with 
experimental observations. They also show that while the 90° spectrum scales with Strouhal number as predicted by 
the Lighthill theory, the small angle spectrum exhibits the experimentally observed Flelmholtz number scaling. 
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Differences in spectral shapes at small and large angles to the jet axis are attributed to the fact that the first 
component of the Green’s function is dominant at 0 = 90° while the second component is dominant at small 0 and 
not to refraction and the spectral broading resulting from a Doppler frequency shift. 

The authors would like to thank Dr. Lennart Hultgren for the use of his Rayleigh equation code and Dr. Louis 
Handler for providing a subroutine used in the spectral computations and Dr. Milo Dahl for his comments on the 
manuscript. 


Appendix A 


As in (5.1) and (5.2), the mean pressure, density, and Reynolds stresses expand as 

p = P(X, x L ) + £ P {1) (X, x ± ) + s 2 R (2) (X, x ± ) + . . . 


(Al) 


and 


P = R(X,x ± ) + eR (1) (X,x ± ) + ... 


(A.2) 


T =T 
w w 


( 0 ) 


+ e 7’(f) +e 2 r 


( 2 ). 


(A.3) 


with similar expansions for h, Hj . and H 0 implied by equations (2.6), and (3.11). Substituting these into the mean 
flow equations (2.5) with cr (/ . = q, = 0 . and assuming that the Reynolds stresses vanish in the free stream, shows that 
the result will only balance if 
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and 
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(no sum on y ^ / = 2, 3) . 


The lowest order mean flow equations become 
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where the operator DJDt is now given by 
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(A. 12) 
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Equations (2.22), (A.l) and (A.3) imply that 
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and it therefore follows from (A.5) and (A. 6) that 
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Appendix B 


The first order perturbation is determined by 
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where the right hand sides are explicitly given by 
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and S is given by (A12). 


Appendix C 

Eliminating the velocity-like variable from (B.l) to (B.3)( as was done for the 0 th order solution) to obtain 
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This result, together with (5.16), shows that g^'J must have a two term decomposition similar to that of (5.21). The 
first term remains uniformly small compared to the corresponding term in the zeroth order solution and need not be 
considered further. The second must be of the form 
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where cp^ (jcj_) is determined from 
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and we have assumed, for simplicity, that the higher order terms in the expansions for the empirically determined 
stresses f ■ have been chosen so that t/^ = c 2 = 0 in order to make independent of both X\ and t. We have set 
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and comparing (C.l) to (C.4) with (B2), (B3), and (5.7) shows that 
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with j = 2,3 and no sum on n. 
Since C k is self adjoint and 


4, 9„/ =0 


(C.8) 


by definition, multiplying (C.3) by (p„/, integrating over x±_ and using the results in appendix A shows that the left 
side of this result must vanish and it follows that A„/ must satisfy (6.3) with 
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where we have used (C.8) and integrated by parts to obtain this result and 
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with j = 2,3 and no sum on /?,/.. 

Equations (B.2) and (B.3) show that (C.10) can be written as 
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where DJDt is now given by (A. 1 1 ), and cp lfrf =-cp ln/ /( y-l) and <p 5nl = (p 5n/ /(y-l) can be obtained from (C. 5) 

and (C.6). 


Appendix D 


To evaluate the second term, note that a must expand like 


a = ^-a 0 (£ s )/r l -a x (Q/r]_ +... asR L — >oo 
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where 
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when 0 and the Mach number Mare sufficiently large. Inserting this into (6.19) and (6.20) and expanding the results 
shows that 


k = K(a 0 + a l /R 1 +...) = K(a c )+ K , (a 0 )a 1 /i ? 1 +... . (D.3) 


Then inserting this together with (D.l) into (6.21) and expanding that result shows that 
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This means that 
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(see (7.4), (7.5), and (D.2)) and 
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Inserting these into (D.3) shows that k behaves like 
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and inserting this into (6.18) shows that when 0 and M are sufficiently large 
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(D. 10) 
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. were R is the radial coordinate defined in (7.6) and 0^ (9,X') is defined by (7.10). 
Equations (6.21), (D.l), and (D.7) imply that 
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and therefore that the amplitude factor in (6.30) behaves like 
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Differentiating (7.10) and using (6.19) to (6.21), (D.2), (D.3), and (D.12) shows that 
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which means that 
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Appendix E 


Taking the Laplace-Fourier transform of (5.8) yields 
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where 0 a is defined by equation (2.15) and £ /; denotes the Rayleigh operator (5.11). Transferring the derivatives of 
the 5-function to the source variable x' yields 


( ku -<°) c &o) _ 

□ 2 ^ U ln “ 
c 0 


P2 


(*i) « [kU(x' ± )-cof 8(x 1 -x' 1 ) 


(y -l) [«/(*! )-»]" ” ^(*i) (2n) 


+1 


k8 Gl {kU- oo) R 

i n-> ° a4 

y - 1 cl 




(E.2) 


(2") 


where 5^ is given by equation (2.15) but with x, replaced by x\ and the dependence on the slow variable A has 
been suppressed, since it enters only parametrically. It is now clear that the solutions to these equations can all be 
expressed in terms of the solution G a [co,k, X; x ± of (5. 10) by 
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Figure Captions 


Figure 1 - Turbulent Jet Flowfield 
Figure 2 - Laplace inversion contour. 

Figure 3 - Fourier inversion contour. 

Figure 4 - Plot of -3m0 VS- 0 M =1.1 (solid), M =1.3 (dashed), 

M = 1 .5 (dot-dash); ( a ) n=0 ,(b)n=l,(c)n = 2. 

Figure 5 - Measurements of the far- Held acoustic spectra at 6 = 90 ° and 
0 — 30 f rom Norum (1994). 

Figure 6 - Far-field acoustic spectrum computed from equation (9.4). M = 1.5 ; 

0 = 30 (solid), 0 = 90 (small dashed), expected effect of cross coupling (large 
dashed). 

Figure 7 - Far-field acoustic spectrum. M = 1.5 ; 0 = 30° (solid), 0=40° 
(dashed), 0=60° (dash-dot), 0 = 75° (dotted), $ 0=90° (dash-dot-dot). 

Figure 8 - Far-field acoustic spectrum. M = 1 .3 ; 0 = 30° (solid), 0 = 50° 
(dashed), 0=60° (dash-dot), 0 = 75° (dotted), $ 0=90° (dash-dot-dot). 

Figure 9 - Far - field acoustic spectrum vs. Strouhal number, M = 1.5 (solid), 

M = 1.3 (dash), M = 1.1 (dash-dot) ; ( a ) 0=30° , ( b ) 0 = 90° . 


Figure 10 - Far - field acoustic spectrum vs. Helmholtz number, 0 — 30° ; 
M = 1.5 (solid), M = 1.3 (dash) , M = 1.1 (dash-dot). 
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